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EARLY UNIVERSE 

On the origin of near-infrared 
extragalactic background 
light anisotropy 

Michael Zemcov, 1 ' 2 Joseph Smidt, 3 ' 4 Toshiaki Arai, 5 ' 6 James Bock, 1 ' 2 * Asantha Cooray,* 
Yan Gong, 4 Min Gyu Kim, 7 Phillip Korngut, 2 ' 1 Anson Lam, 8 ' 1 Dae Hee Lee, 9 
Toshio Matsumoto, 5 ' 10 Shuji Matsuura, 5 Uk Won Nam, 9 Gael Roudier, 2 
Kohji Tsumura, 11 Takehiko Wada 5 

Extragalactic background light (EBL) anisotropy traces variations in the total production 
of photons over cosmic history and may contain faint, extended components missed in 
galaxy point-source surveys. Infrared EBL fluctuations have been attributed to primordial 
galaxies and black holes at the epoch of reionization (EOR) or, alternately, intrahalo 
light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report 
new EBL anisotropy measurements from a specialized sounding rocket experiment at 
1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known 
galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely 
explained by IHL emission. The measured fluctuations are associated with an EBL intensity 
that is comparable to the background from known galaxies measured through number 
counts and therefore a substantial contribution to the energy contained in photons in 
the cosmos. 



It near-infrared wavelengths, where the large 
zodiacal light foreground complicates ab- 
solute photometry measurements, the extra- 
galactic background light (EBL) may be best 
accessed by anisotropy measurements. On 
large angular scales, fluctuations are produced 
by the clustering of galaxies, which is driven by 
the underlying distribution of dark matter. EBL 
anisotropy measurements can probe emission from 
epoch of reionization (EOR) galaxies (1-3) and 
direct-collapse black holes (4) that formed dur- 
ing the EOR before the universe was fully ionized 
by exploiting the distinctive Lyman cutoff feature 
in the rest-frame ultraviolet (UV), thus probing 
the UV luminosity density at high redshifts (5). 
However, large-scale fluctuations may also arise 
from the intrahalo light (IHL) created by stars 
stripped from their parent galaxies during tidal 



interactions (6) at redshift z < 3. A multiwave- 
length fluctuation analysis can distinguish among 
these scenarios and constrain the EOR star for- 
mation rate. 

A search for such background components 
must carefully account for fluctuations produced 
by known galaxy populations. Linear galaxy clus- 
tering is an important contribution to fluctuations 
on scales much larger than galaxies themselves. 
On fine scales, the variation in the number of 
galaxies produces predominantly Poissonian fluc- 
tuations, with an amplitude that depends on the 
luminosity distribution. Anisotropy measurements 
suppress foreground galaxy fluctuations by mask- 
ing known galaxies from an external catalog. 

The first detections of infrared fluctuations in 
excess of the contribution from known galaxies 
with the Spitzer Space Telescope (7-9) were in- 



terpreted as arising from a population of faint 
first-light galaxies at z > 7. The Hubble Space 
Telescope was used at shorter wavelengths (10) 
to carry out a fluctuation study in a small deep 
field but did not report fluctuations in excess of 
known galaxy populations. Measurements with 
the AKARI satellite (11) show excess fluctuations 
with a blue spectrum rapidly rising from 4.1 |j.m 
to 2.4 |j.m. Fluctuation measurements in a large 
survey field (6) with Spitzer agreed with earlier 
measurements (7-9) but were instead interpreted 
as arising from tidally stripped stars at z ~ 1 to 3. 
Most recently, a partial correlation has been re- 
ported (12) between Spitzer and soft x-ray im- 
ages, which Yue et ol. (4) interpret as arising from 
direct-collapse black holes at z > 12. 

We have developed and flown the specialized 
Cosmic Infrared Background Experiment [CIBER 
(13)], a rocket-borne instrument specifically de- 
signed to study the spatial and spectral properties 
of the EBL. The imaging instrument (14) mea- 
sures fluctuations in Ak/X - 0.5 bands centered 
at 1.1 and 1.6 urn using two 11-cm telescopes each 
with a 2° by 2° field of view. Here, we report our 
analysis of data from two flights in 2010 and 2012. 

The CIBER imager data are reduced from raw 
telemetered time streams, flat-field-corrected based 
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on a laboratory measurement, and masked for 
stars and galaxies using the Two Micron All Sky 
Survey (2MASS) J- and #-band catalogs (15). We 
analyze differences between fields to reduce the 
effect of flat-fielding errors. The auto- and cross- 
power spectra of the masked, differenced images, 
corrected for mode coupling from the mask using 
a correction matrix, are shown in Fig. 1. We also 
compute auto- and cross-spectra from a Spitzer 
infrared array camera (IRAQ 3.6-|im image that 
coincides with two of the five CIBER fields. 

The CIBER auto- and cross-spectra show a sig- 
nificant excess over the predicted fluctuations 
from known galaxy populations (16) at I < 5000, 



where the multipole moment I = n/ 6, and 6 is the 
angular separation of two points on the sky. The 
excess is also evident in the cross-power spectra 
with Spitzer, showing that the source of the fluc- 
tuations is largely common from 1.1 to 3.6 urn. 
The large-scale fluctuations measured with CIBER 
correlate between the two flights in all combina- 
tions of bands and are independent of the detector 
arrays; the data pass multiple internal consist- 
ency tests (see the supplementary materials). 

We rule out the following sources for pro- 
ducing the large-scale fluctuations (Fig. 1): (i) sun- 
light scattered by interplanetary dust, (ii) starlight 
scattered by interstellar dust (known as diffuse 



Galactic light), and (ii) fluctuations from faint 
stars. Zodiacal light fluctuations are eliminated 
because we observed the same fields separated 
by 17 months and obtained consistent results, 
viewed through two independent lines of sight 
through the interplanetary dust cloud. Constraints 
on zodiacal light fluctuations from Spitzer at 
8 |j.m (8), scaled by the zodiacal spectrum to near- 
infrared wavelengths (17), lie significantly below 
the detected fluctuations (Fig. 1). The zodiacal 
spectrum is scaled by a conservative upper limit 
to the fluctuation amplitude measured at 7 [im 
from AKARI (17, 18) in Fig. 2. 

We estimate diffuse Galactic light fluctuations 
(Figs. 1 and 2) by scaling the Infrared Astronomical 
Satellite (IRAS) 100-|im intensity by a factor con- 
sistent with previous diffuse Galactic light measure- 
ments (19). We estimate the diffuse Galactic light 
component by calculating the CIBER-IRAS 
cross-spectra and fitting a single amplitude coef- 
ficient to a Ci oc Z~ 3 spectrum with a fixed 1.1 urn / 
1.6 |j.m color ratio, which is less than the detected 
CIBER power spectra at all spatial scales. Fluc- 
tuations from the extinction of the background 
light by galactic dust are also negligible. 

We constrain fluctuations from unmasked stars 
using the UKIDSS-UDS stellar catalog (20), which 
is complete to J = 24.9, H = 24.2 (5 o), accounting 
for > 99.9% of the integrated light from stars. We 
compute the autopower spectrum of residual 
stars below our cutoff flux and find that the 
amplitude is negligible on large scales and fol- 
lows a Poisson spectrum. 

The root mean square (RMS) fluctuation am- 
plitude 8AJx = + l)Ci/2n)V 2 over 500 < I < 
2000 has a spectral energy distribution that is 
approximately Rayleigh- Jeans (Fig. 2). The 1.1-um 
data point lies 2o below the best-fitting X~ 3 
scaling, suggesting the possibility of a departure 
from the Rayleigh- Jeans spectrum at short wave- 
lengths. To fit models to CIBER and Spitzer data, 
we mask Spitzer data to a depth L < 16 that is 
comparable to CIBER. The Spitzer power spectra 
give a non-Poissonian I 17 signal at high multi- 
poles for the L < 16 flux cut (Fig. 1), which is 
evidence for nonlinear clustering that is not ob- 
served with a deeper flux cut (6) and is not pre- 
dicted in linear galaxy clustering (16). 

Having excluded explanations based on zodiacal 
light, diffuse Galactic light, and the clustering of 
known galaxies, we consider three possibilities 
to explain these measurements: EOR galaxies, 
EOR black holes, and IHL. For first-light galaxies, 
we use models for population II and population 
III stars (21) and combinations of the two. The 
direct-collapse black hole model at z > 12 is ex- 
cluded both by amplitude and color, as it has no 
mechanism to generate substantial 1.1- and 1.6-|am 
power (4). For IHL, we use models based on the 
spectral energy distributions of old stellar pop- 
ulations in dark matter halos (6). We simulta- 
neously fit the 1.1-, 1.6-, and 3.6-|j.m autopower 
spectra, taking into account galaxy clustering, dif- 
fuse Galactic light, and residual flat-field errors 
at low I. The IHL model, summarized in Fig. 1, 
generally fits the data but somewhat underpre- 
dicts the amplitude at 1.1 and 1.6 |im; neither the 
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Fig. 1. CIBER and Spitzer auto- and cross-spectra at 1.1, 1.6 and 3.6 |im. We show the CIBER auto- 
spectra for 1.1 by 1.1 |am, and 1.6 by 1.6 |am (A), the CIBER 1.1 by 1.6 urn cross-spectrum (B), the CIBER- 
Spitzer 1.1 by 3.6 |im and 1.6 by 3.6 cross-spectra (C), and the Spitzer 3.6 by 3.6 urn auto-spectra (6) 
(D). At 1.1, 1.6, and 3.6 |im, we indicate previous measurements [open circles; (6, JO); note that the Hubble 
measurements apply a much deeper flux cut at 1.1 and 1.6 urn for masking and that the Spitzer flux cut is 
somewhat deeper than the cut we are applying at 3.6 urn]. The 3.6-um points use the same data set but are 
masked to a lower source flux for comparison to our spectrum, which is masked to L < 16. The increased 
depth reduces some of the mid-Z power. We show constraints on astrophysical foregrounds, including 
unmasked stars and z < 5 galaxies (16), zodiacal light (8), and diffuse Galactic light. In all cases, we 
detect a significant excess power at Z < 5000 (angular separations 9 > 4.3')- We model the data 
using components from IHL (6) and z > 7 first galaxies (3, 5), obtaining a 95% confidence upper 
limit on the EOR contribution. The fitted total indicated by the filled band includes all of the astro- 
physical components plus a bounded systematic error for flat-field variations. The width of the band 
indicates the 68% uncertainty interval of the fit, including all of the modeling uncertainties. The cor- 
relation coefficients for 1.1 by 1.6 ^m, 1.1 by 3.6 |xm, and 1.6 by 3.6 are 0.76 + 0.10, 0.55 ± 0.14 
and 0.31 ± 0.14, respectively, with no statistically significant angular dependence in the correlations 
(see section 6.2 of the supplementary materials for a description of this calculation). 
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Fig. 2. The electromagnetic spectrum of the near- 
infrared fluctuations. We show measurements of the 
fluctuation power from CIBER and Spitzer averaged 
between 500 < Z < 2000 (solid points). Also indicated 
are previous measurements from AKARI (11) and 
Spitzer (6, 7) at Z = 3000 that use deeper masking 
thresholds. In all cases, we subtract the contribution 
from the shot noise of unmasked galaxies (16). We 
indicate the best-fitting Rayleigh-Jeans spectrum to 
the points from this analysis, estimates for diffuse 
Galactic light fluctuations (19), a conservative con- 
straint on zodiacal light fluctuations (17, 18), and an 
upper limit on Galactic emission (27). The known 
foreground components have both smaller ampli- 
tudes and different spectra than the measurements. 
We show the residual from the best-fitting Rayleigh- 
Jeans spectrum 8 RJ in the bottom panel, scaled by A 3 
to reduce the range. The shortest wavelength mea- 
surement at 1.1 |j.m is 2.0 o below the fit, indicating a 
possible short-wavelength departure from a Rayleigh- 
Jeans spectrum. 
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Table 1. Contributions to near-infrared EBL anisotropy and intensity. At each wavelength, we list the measured fluctuation amplitude at large angular 
scales; the model-dependent ratio of EBL intensity to EBL anisotropy; the IGL determined by previous measurements; the ratio of the IHL and IGL intensities; 
and finally, the inferred total background intensity from both components. We also list the background intensity that would arise assuming the measured 
fluctuations are entirely due to high-redshift EOR galaxies. 

W..m^ Measured 6kl x * XJ^ihl Mx,iHi_t Mx,igl§ >^x,ihl Mx,\hl + >A,igl Mx,eor\\ 



Mumj (nWm- 2 sr _1 ) 6U X (nW m -2 sr _1 ) (nW m -2 sr _1 ) X/xjgl (nW m -2 sr _1 ) (nW m -2 sr _1 ) 

LI L4lp 5 TO^I 97tf? 07 28 

1.6 1.9±8s 6 11.4if3 9.0t! 2 7 6 1.3 20.4 + 5 6 i° 38 

2.4 0.32±0.05t 7 2.2 + 0.4 7& + -ih 0.3 10.0t 2 3° 6.4 
3.6 0.072too2i 9 °- 65 -Qi| 5.2 ± 1.0 0.1 5.9 + 1.0 1.4 
3.6# 0.049 + o 0 b°07 9 0.44+8o6 5.2 + 1.0 0.1 5.6 + 1.0 1.0 

4.5 0.053 + 0.0231 7 0.37 + 0.16 3.9 + 0.8 0.1 4.3 + 0.8 1.0 



*RMS fluctuation amplitude computed as averages of measured data over 500 < I < 2000, except for those marked f, which are determined at I = 3000 using 
fainter mask cuts due to restricted field size (see also note "). +The IHL background from the product of columns 2 and 3. §The IGL background as 
compiled by (28). ||Computed EOR background assuming EOR fluctuations with = 20. ^Determined at K band corresponding to 2.2 ^m. #Computed 

using the measurements of (6) averaged over 500 < I < 5000. 



IHL nor the EOR models quite match the high 
observed 1.6-/3.6-|jm color ratio (see section 10 of 
the supplementary materials). This may indicate 
that additional components are reflected in the 
data. However, more theoretical work is required 
to determine whether adding nonlinear galaxy 
clustering and nonlinear IHL production to the 
model can improve the fit. 

We estimate the EBL intensity associated with 
the fluctuations by taking the measured fluctua- 
tion amplitude between 500 < I < 2000, obtained 
by subtracting estimated contributions from \ow-z 
galaxies and diffuse Galactic light, and multiply- 
ing by a model-dependent contrast factor }J X / 



bUx, where U x is the total intensity associated 
with a component. For the IHL model, which has 
a low contrast factor, we obtain an associated EBL 
of 7.0+ 4 5 ° and 11.4+ 4 54 nW m -2 sr _1 at 1.1 and 
1.6 jam, respectively. As shown in Table 1, the IHL 
background is of a similar magnitude to the 
integrated galaxy light (IGL) background derived 
from galaxy counts. However, we note that the 
IHL background has a much bluer color than the 
IGL background. We similarly estimate the IHL 
background at longer wavelengths from AKARI 
and Spitzer. Nonlinear galaxy clustering ap- 
pears to contribute to the Spitzer fluctuations, so 
we quote two values that depend on the choice of 



flux cut, the deeper flux cut being less sensitive 
to the nonlinear clustering contribution. The 
CIBER data do not appear to be as sensitive to 
the flux cut, perhaps due to the higher IHL to 
IGL ratio at these wavelengths. 

The total EBL, the summation of the IHL and 
IGL backgrounds, is consistent with current EBL 
measurements. Near-infrared absolute photometric 
background measurements remain uncertain due 
to the bright zodiacal foreground, but the lowest 
such measurement (22) gives 21 ±15 nW m" 2 sr" 1 
and 13.3 + 2.8 nW m" 2 sr" 1 at 1.25 and 3.6 |^m, 
respectively. The High Energy Stereoscopic Sys- 
tem (HESS) measurement of the EBL using y-ray 
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absorption spectra (23), 15 + 2 (statistical) + 3 
(systematic) nW m" 2 sr" 1 at 1.4 jam, is indepen- 
dently consistent with the sum of the IHL and 
IGL backgrounds, within the uncertainties in 
all measurements. 

An EOR interpretation of these measurements 
would result in a large background, as the contrast 
factor is larger than that of IHL, 10 to 100 for EOR 
models (3, 6). The implied EOR background (see 
Table 1) originating from high redshifts is diffi- 
cult to justify due to the overproduction of me- 
tals and x-ray background photons (24). 

Our results indicate that a substantial fraction 
of the EBL at optical and near-infrared wave- 
lengths originates from stars outside of galaxies 
(with boundaries as traditionally defined). This 
in turn adds to the cosmic energy budget and, 
depending on the mass characteristics and spec- 
trum of the population responsible, could help 
alleviate the "photon underproduction crisis" 
(25) and the "missing baryon problem" (26). We 
see no evidence for a detected EOR background 
component in our data. Multiwavelength fluctu- 
ation measurements extending into the optical 
will help discriminate the EOR background com- 
ponent using the redshifted Lyman cutoff, and 
future spectroscopic measurements will enable 
tomographic measurements to determine the 
history of IHL production. 
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EARLY EARTH 

Sulfate was a trace constituent of 
Archean seawater 

Sean A. Crowe, 1 ' 2 *! Guillaume Paris, 3 * Sergei Katsev,* CarriAyne Jones, 1 ' 2 

Sang-Tae Kim, 5 Aubrey L. Zerkle, 6 Sulung Nomosatryo, 7 David A Fowle, 8 Jess F. Adkins, 3 

Alex L. Sessions, 3 James Farquhar, 9 Donald E. Canfield 2 

In the low-oxygen Archean world (>2400 million years ago), seawater sulfate 
concentrations were much lower than today, yet open questions frustrate the translation 
of modern measurements of sulfur isotope fractionations into estimates of Archean 
seawater sulfate concentrations. In the water column of Lake Matano, Indonesia, a 
low-sulfate analog for the Archean ocean, we find large (>20 per mil) sulfur isotope 
fractionations between sulfate and sulfide, but the underlying sediment sulfides preserve 
a muted range of 5 34 S values. Using models informed by sulfur cycling in Lake Matano, 
we infer Archean seawater sulfate concentrations of less than 2.5 micromolar. At these 
low concentrations, marine sulfate residence times were likely 10 3 to 10 4 years, and 
sulfate scarcity would have shaped early global biogeochemical cycles, possibly restricting 
biological productivity in Archean oceans. 



Sulfur interacts with carbon and oxygen in 
global biogeochemical cycles that regulate 
Earth's surface chemistry and biology (1). 
At 28 mM, sulfate is abundant in modern 
seawater, fueling extensive sedimentary mi- 
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crobial sulfate reduction (MSR) (2). At these con- 
centrations, MSR typically imparts large sulfur 
isotope fractionations (3), allowing the use of 
sulfur isotopes to reconstruct past global change 
(4, 5). Small sulfur isotope fractionations preserved 
in bulk pyrite from Archean (>2400 million years 
ago) rocks led to the original conclusion that 
the Archean oceans contained <200 uM sulfate, 
or ~1% of modern seawater (5). The distribution 
of mass-independent sulfur isotopes in Archean 
sediments (4, 6-8) and box models of global sulfur 
cycling (9) imply even lower Archean seawater 
sulfate of <60 to 80 uM. Paradoxically, microscale 
sulfur isotope data from Archean pyrites (10-12) 
reveal large sulfur isotope fractionations of up to 
40 per mil (%o) in the Archean— fractionations 
only seen at hundreds to thousands of micro- 
molar sulfate in modern environments (13- IS). 
Electron-donor availability (3, 16) and MSR rate 
(3, 16, 17), however, also exert influence, with larger 
fractionation typically imparted when electron do- 
nors limit sulfate-reduction rates, allowing the 
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